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( An TV-channel generalization of the network model of Chalker and Coddington is considered. 

The model for TV = 1 is known to describe the critical behavior at the plateau transition in systems 
exhibiting the integer quantum Hall effect. Using a recently discovered equality of integrals, the 
f"S i network model is transformed into a lattice field theory defined over Efetov's a model space with 

• unitary symmetry. The transformation is exact for all TV, no saddle-point approximation is made, 

and no massive modes have to be eliminated. The naive continuum limit of the lattice theory is shown 
to be a supersymmetric version of Pruisken's nonlinear a model with couplings a xx — a xy — TV/2 
at the symmetric point, ft follows that the model for TV = 2, which describes a spin degenerate 
Landau level and the random flux problem, is noncritical. On the basis of symmetry considerations 
and inspection of the Hamiltonian limit, a modified network model is formulated, which still lies 
in the quantum Hall universality class. The prospects for deformation to a Yang-Baxter integrable 
, vertex model are briefly discussed. 
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I. INTRODUCTION 



The single-electron states of a two-dimensional disordered electron gas in a strong magnetic field are localized 
except at the energies of the Landau band centers. As the Fermi energy approaches such a band center, a critical 
phenomenon takes place: the localization length diverges and the Hall conductance jumps from one plateau to the 
next. This phase transition, which belongs to the general class of Anderson metal-insulator transitions, has been 
(f) , seen in several experiments and a substantial amount of data on its critical behavior is available from a number of 
^ ■ numerical simulations (in the absence of electron-electron interactions), see |jj and references therein. Unfortunately, 
in spite of considerable efforts expended over the last decade, our analytical understanding of the plateau-to-plateau 
transition is still rather poor. It is expected that the critical behavior is described by some nonunitary conformal field 
' theory, but this theory has not yet been identified. 

There exist two opposite limits from which the transition in the noninteracting system can be approached theoret- 
' ically. The first limit is that of a slowly varying random potential with a correlation length l c much larger than the 
^\ , magnetic length Is- In this limit, the electron's motion can be described in semiclassical terms More precisely, the 
motion separates into a rapid cyclotron motion superimposed on a slow guiding center drift along spatially localized 

Sequipotential lines. As the Fermi energy approaches the center of a Landau band, a percolating path develops and a 
localization-delocalization transition takes place. Close to the transition, the quantum mechanical possibility for an 
electron to tunnel from an equipotential to a neighboring one is a relevant perturbation. The essential features of 
this quantum percolation transition were cast into a random network model by Chalker and Coddington || . In their 
model an electron acquires random U(l) phases while moving along the directed bonds of a square network, and is 
£j , scattered to the right or left every time it passes through a node of the network. The model will be reviewed in more 
IV ' detail below. Suffice it to say here that the model has been studied by numerical simulation but has, in its original, 
. £h ! spatially isotropic formulation, defied analytical treatment up to now, although a certain amount of analytical insight 
■ has come from consideration of its anisotropic j| and weak disorder || versions. 

The opposite limit is l c <C Ib- Historically, this was the first limit to be studied analytically, the reason being 
that it is this limit that was amenable to the field theoretic machinery developed at the beginning of the 1980's by 
Wegner ||, Efetov |?J and others. Starting from a Gaussian white-noise potential (l c = 0) Pruiskcn ||] used the 
replica trick to set up a generating functional for the disorder averaged retarded and advanced Green's functions of 
the single-electron Hamiltonian. After a Hubbard-Stratonovitch transformation to matrix-valued fields, he made a 
saddle-point approximation, valid in the limit of a high Landau level. This was followed by a gradient expansion 
leading to a U(n+ + n_)/U(n + ) x U(n_) nonlinear a model with a vanishing number of replicas, n+ = n_ = 0, and 
a parity-violating topological term due to the strong magnetic field. The coupling constants of the model, a xx and 
a xy , were identified Q with the physical conductivities of the 2d disordered electron gas. The topological coupling 
9 = 2ira xy , by its very nature, has no effect on the equations of motion of the classical field theory, but does change 
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both the Hamiltonian and the symplectic structure and, consequently, the commutator of the quantum theory. It was 
argued fic| that the nonlinear a model, while generically being massive (i.e. having a finite correlation length) in two 
dimensions, becomes critical at — it. The vanishing of the mass gap corresponds to the appearance of delocalized 
states at the center of a Landau band. Thus, Pruisken's model provided the right kind of scenario in which to develop 
a scaling theory of the plateau-to-plateau transition. A supersymmetric version of the model first appeared in JlT| . 

Undeniably, Pruisken's nonlinear a model has served as a great inspiration to theory. One of its early successes 
was the a xx -a xy flow that was conjectured from it ]l2] , [l0| and later verified by numerical and real experiments, see 
jl| and references therein. In spite of this, Pruisken's model or, rather, its promoters have been criticized, see (jlf l 
for a summary. For one thing, the model has never yielded any quantitative results for the critical behavior at the 
transition, and much less has it been solved. (The same statement applies to the general class of nonlinear a models 
with a topological term. None of these has ever been solved, at least not directly.) For another, even the derivation 
of the model is vulnerable to criticism: the validity of the saddle-point approximation that is made to eliminate the 
so-called massive modes, requires a xx ^S> 1. Although this inequality is satisfied for the bare (or SCBA) value of a xx 
in the limit of a high Landau level, the renormalized theory at a xy — 1/2 is expected to have a a xx of order unity or 
less. The cure proposed by Pruiskcn was to assume the renormalizability of his model, and appeal to the RG flow to 
take the coupling constant a xx from large to small values. However, such an assumption needs to be justified and, 
in fact, is not acceptable by current field theoretic knowledge, for Pruisken's model apparently lacks the conformal 
structure that is required of a fixed point theory with a continuous symmetry in two dimensions. (In other words, 
the model, while definitely being critical at a xy = 1/2, does not possess the conservation laws expected of an infrared 
stable fixed point.) 

Two advances will be made in the present paper. The first is to establish a very close connection between Pruisken's 
nonlinear a model and the network model of Chalker and Coddington at criticality. We will show that the latter can 
be viewed as a lattice discretization of the former or, conversely, taking the continuum limit of the network model 
yields the a model. A more detailed outline is the following. We start out by reviewing the supersymmetric version 
of Pruisken's model and the Chalker-Coddington network model in Sees. || and III. Then, in Sec. IV, the network 



model is reformulated as a lattice-regularized field theory defined over Efetov's supersymmetric nonlinear a model 
space with unitary symmetry. This reformulation is exact. In contrast with the conventional method due to Wegner, 
Efetov and others, no saddle-point approximation is needed to eliminate the massive modes. Moreover, Sec. [v] shows 
that not only is the supersymmetric reformulation of the network model defined over the same field space, but it also 
has the same global symmetries as Pruisken's model. At the critical point of the network model, where the correlation 
(or localization) length diverges, the long wave length physics of the supersymmetric lattice theory is expected to be 
described by a continuum field theory. The symmetries dictate that this continuum theory be Pruisken's model at the 



critical coupling a xy — 1/2. As is arg ued in Sec. VIIA, the other coupling constant, a xx , equals 1/2. The numerical 
value of <T xy is checked in Sec. VII B, by evaluating the lattice action on a smooth field configuration with nonzero 
topological charge. Sec. VII C extends these results to an Af-channel network model with random U(N) matrices on 
links. 

What we learn from all this is that, although historically Pruisken's model was first obtained from the Gaussian 
white noise limit l c <C Ib, in a low Landau level it is actually more closely related to the opposite limit l c ^> Ib, since 
it is the latter that provides the microscopic basis for the network model. 

The mathematical basis underlying the exact transformation from the A-channel network model to the super- 
symmetric lattice field theory is quite natural and simple, and is briefly sketched as follows. (Readers who are not 
interested in mathematical structures may want to skip this paragraph.) For a pair of positive integers n, N consider 
the tensor product C" ® C , on which the group GL(nN) acts by linear transformations. Span the corresponding Lie 
algebra gl(niV) by bilinears {ca^bYabI^i' n m fermionic creators c and annihilators c, which act in a Fock space with 
vacuum |0). There exist two natural subalgebras: gl(n), generated by J2i^A c B^ an< ^ gK-^O; generated by YIa^A^a- 
Now put n = n- + n + and fill the "negative energy states" to form |vac) = ni=i IL4=i The particle-hole 

coherent states [Q that are generated by the action of GL(n) on |vac), are holomorphic sections of a line bundle 
associated to the homogeneous space G/H := GL(n)/GL(n + ) x GL(n_) by the Slater-determinant representation of 
GL(n_) on |vac). They are parameterized by a complex n + x n_ matrix Z with adjoint Z'. By combining the closure 
relation for particle-hole coherent states with a few elementary properties of Fermi-coherent (or Grassmann-coherent) 
states, one can prove the following equality of integrals: 

n+ 



[ dU TT Det (1 - e +itp+a U) f] Det (l - e^- 6 t/ f ) 

Mm a=l 6=1 

J dp(Z, Z f ) Bet- N (1 + Z^Z) T)ct N (l + tfe^+Ze*- 
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where ip± = diag(y?±i, <p± n ±) are diagonal matrices with real entries, dU is the Haar measure of a U(iV) subgroup 
of GL(iV), and d[i[Z 1 Z') expresses the U(n)-invariant measure of a compact symmetric space U(n)/U(n + ) x U(n_) 
contained in the complex homogeneous space G/H. This integral identity forms the mathematical basis of our 
formalism. Its supersymmetric extension ]l5| permits to carry out the disorder average over the random U(iV) 
matrices placed on the links of the network, at the expense of introducing an integration over fields Z taking values 
in a symmetric superspace. 

Another issue addressed in this paper is the question whether the supersymmetric formulation of the network model 



offers the possibility for an exact analytical solution. Sec. VI reveals that the model resembles what is called a vertex 



model in statistical physics, in the sense that the Boltzmann weight is a product of factors, one for each node, or vertex. 
The w eight associated with a single vertex is called the R-matrix. The symmetries of the R-matrix are investigated 



in Sec. VIII. It is eventually found that it can be interpreted as a map TZ:V<Si)V—>V®V, where V is an irreducible 
lowest-weight module for the Lie superalgebra gl(«, n), and n = n + + n_ = 1 + 1 for the case of one retarded and one 
advanced Green's function. This looks interesting as one may hope that 1Z can be deformed to an R-matrix that solves 
the quantum Yang-Baxter equation underlying the integrability of two-dimensional vertex models. One would then 
have the possibility of an analytical and exact computation of critical properties. Unfortunately, the specific choice 
of local directions for the single-particle motion on the network, shown in Fig. 2 below, turns out to be incompatible 
with the standard schemes jl6| for constructing solutions of the quantum Yang-Baxter equation. The reason is that 
the R-matrix of a Yang-Baxter solvable model always transfers from one side of the vertex to the other, whereas the 
Chalker-Coddington vertex maps the horizontal degrees of freedom into the vertical ones, or vice versa. Thus the 
Chalker-Coddington model in its original form does not fit into the canonical framework of the theory of integrable 
systems, and I am not aware of any method to make analytical progress with it. 

However, this is not yet the end of the story. Additional insight can be gained by considering the anisotropic limit 



m of the Chalker-Coddington model. This limit and its relation to Pruisken's model are reviewed in Sec. [EX. Based 
on it, Sec. [X] proposes a modified version of the isotropic one-channel network model, which differs from the original 
one in two respects. First, the direction of the single-particle motion does not alternate constantly between being 
horizontal and vertical. Instead, an electron may either pass straight through a node (with no change of direction), or 
else be scattered either to the right or to the left. Analysing such a model by the mapping onto Pruisken's nonlinear a 
model, we find that it is likely not to be in the quantum Hall universality class, but in a massive (Haldane type) phase. 
Therefore, a second modification is proposed, which is to add a second channel of propagation on half the links, say 
the horizontal ones, see Fig. 5b below. By the mapping onto Pruisken's model, the doubly modified network model is 
expected to be critical in a range of values of the parameter characterizing the scattering at the nodes. Moreover, by 
the changed transfer dynamics, the standard schemes for constructing solutions of the quantum Yang-Baxter equation 
are no longer ruled out. I hope to elaborate on this theme in a future publication. 



II. SUPERSYMMETRIC FORMULATION OF PRUISKEN'S MODEL (DEFINITIONS) 

The original formulation ^] of Pruisken's nonlinear a model relied on the use of the replica trick. When applied 
to phenomena that are nonperturbative in the disorder strength, the replica trick is not mathematically sound but 
has been demonstrated to lead to incorrect results, at least in some instances |17|] . (The analytic continuation to 
a vanishing number of replicas is not unique in general.) Fortunately, we can avoid the replica trick by using an 
alternative, supersymmetric formalism which is on firm mathematical ground. The purpose of this section is to 
briefly review the supersymmetric formulation fll|| of Pruisken's model in a language that is well suited for what will 
follow below. For simplicity, only the model pertaining to one retarded and one advanced Green's function is treated. 
A more detailed discussion of the model can be found in Jll| . 

To define Pruisken's nonlinear a model we first specify its field space, as follows. Consider a pair Z, Z of 2 x 2 
complex supermatrices 

Z = ( Zbb Zbf I Z = I "? BB "? BF I 

V Zfb z-py j ' y Zfb Zpp J ' 

where the subscripts B and F stand for Bosonic and Fermionic, and let g = ( ^ ^ J S GL(2|2) act on these by 



Z^> g ■ Z = {AZ + B)(CZ + D) -1 , 
Z^ g ■ Z = (C + DZ)(A + BZ)- 1 . 
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Following |l3| , one identifies the pair Z, Z as a set of coordinates for the complex coset space G/H :— G1(2|2)/GL(1|1) x 
GL(1|1), where the denominator is the subgroup generated by the block-diagonal matrices h = ( l5 2. ) . With this 



x D / 

identification, the action Z i— > g-Z and Z ^ g-Z coincides |l3| with the natural action of G on G/H by left translation. 
The supermatrices Z and Z then are viewed as left-translates of the origin in G/H, by writing Z = g ■ = BD~ X and 

The coset space G/H, being a homogeneous space, admits only one (up to multiplication by a constant) rank-two 
supersymmctric tensor that is invariant under the action of G — GL(2|2). In the coordinates Z, Z it is given by 

g = STr(l - ZZ)- 1 dZ(l - ZZ)~ x dZ, 

where STr denotes the supertrace. The G-invariant superintegration measure that derives from g is denoted D(Z, Z) 
and has the localf] coordinate expression ]T^ | 

D(Z, Z) = dZ BB A dZ B B A dZ FF A dZ FF 



dZBFdZBFdZpBdZ-pB 

The integration domain for the bosonic variables is fixed by the conditions 

Zff = —Zff, Zbb = +Zbb, and |Zbb| 2 < ^-i 

where the bar denotes complex conjugation. These conditions select a submanifold Mb X Mf of G/H, 

Mb = U(l, 1)/U(1) x U(l) ~ H 2 , Mf = U(2)/U(l) x U(l) ~ S 2 , 

on which the metric g is Riemann. The variables Zff and Z BB can be shown |l3f| to have a meaning as complex 
stereographic coordinates on the two-sphere S 2 and two-hyperboloid H 2 , respectively. The triple (G/H, Mb x Mp,g) 
is a Riemannian symmetric superspace pQ] , and is called Efetov's cr model space with unitary symmetry. In most of 
the existing literature, this space is parameterized by a 4 x 4 supermatrix Q, introduced below. 

After these preparations, Pruisken's nonlinear a model (or, rather, the supersymmetric version thereof) is defined 
by the functional integral 

(.) = Jv(Z,Z) . exp-S cont [Z,Z] , 
where S con t[Z, Z] = J Ld 2 r is obtained by integrating the following Lagrangian: 

C„ v = STr(l - ZZy^d^Zil - ZZ)- x d v Z, (1) 
and the functional integration measure is 

V{Z,Z)= H d(z(v),Z{v) 

r=(x,y) 

By construction, this field theory is invariant under global Gl(2|2) transformations Z(r) i— > g-Z(r) and Z(y) i— * g- Z(r) 
(each of the terms C^ v is). The partition function Z = (1) equals unity by supersymmetry. To obtain physical 
observables, one adds sources and takes functional derivatives, as usual. 

For some purposes it is useful to switch to a coordinate-independent language by introducing a 4 x 4 supermatrix 
field Q by 



Q = 



1 Z\ ( 1 OWl Z" 



Z 10 -1 \ z 1 



1 Superintegration measures, also called integral superforms or Berezin forms, generically suffer from a coordinate ambiguity, 
or anomaly, see [19). A general procedure by which to define D(Z, Z) globally is described in Sec. II A of [20|. 
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FIG. 1. The network model of Chalker and Coddington in a quasi Id geometry. A single electron propagates ballistically 
along the links and is scattered at the nodes of the network. 



in terms of which the Lagrangian takes the familiar form || 

C xx + C yy = STr (d x Qd x Q + d y Qd y Q)/8 =: L , 
£xy — £yx — STr(d x Qd y Q — d y Qd x Q)Q /8 =: Ltop- 

The global action of GL(2|2 ) on Q is Q(v) i— > gQ(r)g^ 1 . On a configuration manifold C without boundary the integral 
of -L top is quantized [ ]10{Jlq ]: J c L top d 2 r — 2-7rm, and the integer n is called the "winding number" or "topological 
charge" . Lq and L top are invariant under rotations in the xy plane, and these are the only invariants that can be 
formed from £ M „. Note that L is real- valued, whereas Ltop is purely imaginary. The coupling constants a x °J and 

a x °J have been interpreted Q as the longitudinal and Hall conductivities of the 2d electron gas. 
Let me mention in passing that the terms in the Lagrangian can be written as 

L °= (aS? + w) insDct - z{x ^ x '^y L w ■ 

where SDet denotes the superdeterminant. The function lnSDct(l — ZZ) is called the Kahler potential of the Kahler 
supermanifold G/H parameterized by the complex coordinates Z and Z. 



X—T/ .ll — ll' 



III. THE CHALKER-CODDINGTON MODEL 

The Chalker-Coddington model was formulated || for the purpose of describing the plateau-to-plateau transition 
in systems exhibiting the integer quantum Hall effect. Its microscopic foundation as a model for the motion of a single 
2d electron subject to a smooth random potential and a strong magnetic field, is explained in the original paper. 

The model is defined on a finite or infinite square lattice forming a network of nodes and directed links, see Fig. 1. 
A "wave function" of the model is a collection of complex amplitudes one for each link I of the network. 

The dynamics is governed not by a Hamiltonian but by a unitary operator U , called the one-step time evolution 
operator. U acts on wave functions, and is given by a sequence of two distinct operations. The first one is stochastic 
and, in a microscopic picture, accounts for the random phase acquired by the electron's guiding center while drifting 
along equipotentials between saddle points of the random potential. In formal language, the stochastic phases are 
encapsulated in a unitary operator Uq that is diagonal in the basis provided by the links: Uq(V , I) = 5(1', I) expz(p(Z), 
where tp(l) are uncorrelated random variables with a uniform distribution on the interval [0, 27r]. 

The second process built into the one-step time evolution operator is deterministic (i.e., nonrandom) and accounts 
for the quantum mechanical possibility for an electron to tunnel across a saddle point. This is modeled via the nodes 
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FIG. 2. Chalker-Coddington network with two interior contacts M and N. The corresponding links Im and Zjv are connected 
to one incoming and one outgoing channel each. 

of the network. One imagines that the guiding center drift motion of the electron follows the direction indicated by 
the arrows in Fig. 1. An electron incident on a node can be scattered either to the left or to the right, corresponding 
to the two possibilities of continuing on its way along the same equipotential, or tunneling to a neighboring one. It 
cannot be backscattered, and it cannot pass straight through a node. The probabilities for scattering to the right or 
left are denoted p and 1 — p. The probability amplitudes for scattering at the nodes then have magnitude ^fp and 
y/1 — p respectively. They are taken to be real, but they cannot all be chosen positive as this would violate unitarity. 
Various choices are possible. For definiteness we take the amplitudes to be negative for right-up to left-up and for 
right-down to left-down turns (see Fig. 1) and positive otherwise. (Which choice is made actually turns out to be 
immaterial. All we need is that some choice consistent with unitarity exists.) All this fixes a unitary operator Ui 
with matrix elements Ui(l',l) that vanish unless the scattering process / — > I' is allowed by the (one-step) dynamics, 
in which case Ui(l',l) = ±^/p or Ui(l',l) — ±^/l — p, as specified. The full one-step time evolution operator is the 
product U — U\Uq. 

For the purpose of doing numerical simulations, one usually takes the network to be a long strip (quasi-ld geometry) 
with periodic boundary conditions for the short direction to minimize finite size effects. By computing the exponential 
growth of the transfer matrix for the strip and averaging over many realizations of the disorder embodied by Uq, one 
extracts a Lyapunov exonent or inverse localization length. For p = 1/2 the localization length is found to grow 
linearly with the width of the strip, indicating a critical point. 

The existence of critical behavior at p = 1/2 can be anticipated by the following argument j|. Consider the 
Chalker-Coddington model with the open boundary conditions of Fig. 1. For p = (left turns only) all electron states 
encircle elementary squares of the network in the counterclockwise direction, and thus are strongly localized. The 
same statement applies for p = 1 (right turns only) , except that now the orientation of the motion around squares is 
clockwise, and an edge state at the boundary of the strip appears, as is seen by inspection of Fig. 1. The appearance 
of an edge state implies that somewhere in the interval p £ [0, 1] a delocalized state must form. For symmetry reasons, 
this is expected to happen at the left-right symmetric point p = 1/2. 

The method to be introduced in Sec. ^ is quite capable of dealing with a quasi-ld geometry and its boundary 
conditions. However, for maximal simplicity we will consider below a somewhat different setup, where periodic 
boundary conditions are imposed in both directions, i.e. the network is placed on a torus S 1 x S . We envisage 
making a conductance measurement between two interior contacts M and N as shown in Fig. 2. These contacts are 
"point" contacts in the sense that they each attach only to a single link, Im and l^. For most purposes it is helpful to 
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imagine that each link with a contact is replaced by two links, an "in-link" feeding the network from a reservoir, and an 
"out-link" draining the network through an outgoing channel leading back to the reservoir. Outgoing-wave boundary 
conditions are imposed on out-links, i.e. probability flux that is scattered into such a link by the action of U%, must 
exit and never returns to the network. The conductance pertaining to two interior contacts M and N can be computed 
from the Landauer-Biittiker formula, gMN — \Smn\ 2 , where Smn is the S-matrix element relating an incoming state 
on the in-link of In to an outgoing wave amplitude on the out-link of Im- This S-matrix element is determined by 
the solution of Schrodinger's equation Uip — e la ?p with incoming-wave boundary conditions at N. The eigenphase 
e la may be gauged away by a global shift of the random U(l) factors on links, Uq(1) i— > e %a Ua(l). It is convenient 
to implement the loss of probability flux through out-links by the modification Uq(1) = expi<p(l) — > Uo(l out ) = for 
every link I with a contact. Then, solving the equation -0 = Uip by iteration, we easily see that the d.c. conductance 
equals 

9mn = \Smn\ = (IaTITUn) 
where T is the operator 

T = U 1 + U.UoU, + UxUoUiUoU! + ... = (1 - U^y 1 Uj, (3) 

Note that owing to UoQffi) = UoQffi) = 0, the eigenvalues of U-JJq lie inside the unit circle (rather than on the 
unit circle) in C, so the inverse (1 — UiIIq)^ 1 is well-defined. Without changing the matrix element |Sa/jv| 2 , we se t 
UoK) = U (In) = 0. 

IV. REFORMULATION BY SUPERSYMMETRY 

In this section we will show how to map the problem of calculating the disorder averaged conductance (gMN) on 
an equivalent lattice field theory. From gMN 

mean conductance (gMN) equals the disorder average of the expression 

,. vl (l -[/+[/,!)-%> ul(l m ,l° A f) (l^Kl-UxUo)- 1 ^) Ui(l n ,$). 



(l° A T\T\l%) , T = (1 - IWo)- 1 ^, and !* = (!- L/Jtft)- 1 ^, the 



This disorder average can be processed by a variant of Efetov's supersymmetry method introduced in |lq| , as follows. 
The method starts by expressing the product of matrix elements as a Gaussian superintegral: 

(i N \(i - ulu^iu^fKi u.Uo)- 1 ^) 

I 

x exp ( - [6(1', I) - U x (l', l)U (l)]i/j +(T (l) 

-$-r(l')[5()',l) - UiUtfUaWty-Alj), 

where ip±o- are the components of a complex superfield I i— > ip(l). The index a = B or F distinguishes between 
Bosonic and Fermionic components, and ± relates to retarded (T) and advanced (T') Green's functions. The bar 
denotes complex conjugation. D(ip,ip) denotes the "flat" superintegration measure, i.e., is given by the product of 
the differentials of the bosonic components times the product of partial derivatives with respect to the fermionic ones. 
The summation convention is used here and throughout the paper. 

For the following step, it is notationally convenient to absorb Ui temporarily by setting ij]+v(l) — i[i+ a (l')Ui(V , I) 
and V>_ T (0 = $-t(1')Ui(1,1'). We will now trade the average over random phases J_\~\i dUn(l) — JY\ l dnp(l)/2n 
for an integral over Efetov's a model space with unitary symmetry, reviewed in Sec. ||. This is done by a kind of 
Hubbard-Stratonovitch transformation, 



II du o(l) exp (j> + <,(l)Uo(l)i>+a(l) + ^_ T (0E/b(0lM0 y 
l 

Jj'l>M (Z(1),Z(1)) exp (^ +a (l)Z aT (l)^ T (l) + ^ r (l)Z ra (l)il>+v(l)) , (4) 
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which is a special case of a more general result proved in j[f|. The 2x2 supermatrix valued fields I h-» Z(l) and 
I i— > Z(7) are lattice discretizations of the continuous fields of Pruisken's model. The superintegration measure, 

D/i(Z, Z) := const x Z) SDet(l - ZZ), 

is normalized by the condition J Dfi(Z, Z) = 1. Because of the boundary condition at the contacts, Uq(Im) = 
Uq(In) = (in- and out-links), the product over links on both sides of (Q) excludes Im and In- An alternative scheme 
of implementing the boundary condition is to let the product run over all links, including the boundary ones, and set 

Z(l M ) = Z(l M ) = Z(l N ) = Z(l N ) = 0. (5) 

We adopt this scheme. By using (Q) to deal with the disorder average of the product of matrix elements, we obtain 

{i%\(i - uMr^ii^ni LWonM; 

l[Dfx(z(i),z(i)) fl[D(m,m) ^MmMirn^MiTW+BiQ 



x exp 



( - i> +a (l)^ +a (l) + $ +(T (l')Ui(l', l)Z aT {l)i;_ T {l) 

-^_ T (l)V-r(0 + ^-r(O^l(i.O^(0^(0 



The final step is to carry out the Gaussian integral over ip and ip, which gives 



(qmn ) = 



f*[[DiJi(z(!),Z(lj) SDet«(l-ZtfZ) 
^ i 

xlZuil-ZZu)- 1 } (llfJl^xlz^il-ZZ^)- 1 } (Q,l%) 

L J BB L J BB 



by an elementary calculation. Here Z\j denotes the supermatrix field Z evolved forward in time by one action of the 
deterministic scattering operator {Ji, i.e. Zjj — U 1 ZLf\ or, in more explicit notation, 



(u.ZUl) (1,1') = U x (l,l")Z GT {l")U x {l\l"). 

\ /err 



Similarly, Z v ^ is the supermatrix field Z evolved backward in time by one step: Zjj\ — \j\ZU-y. The subscript Tt 
on SDet indicates that the superdeterminant runs over both superspace and the Hilbert space of wave functions 
supported on links. 

In summary, we have shown that the average conductance pertaining to two interior contacts M and N can be 
computed as a two-point function (guN) = (O(M)O(N)) of operators 

O(M) = [Z V (1 - ZZu)- 1 ] (Iff, Iff), 
6(N)= IZu^l-ZZut)- 1 ]^,^), 



in a lattice field theory 



[Y[D(Z(!),Z(1)) • exp-SUM, 



with lattice action 



Sutt[Z, Z] = lnSDet«(l -Z V Z)- lnSDet w (l - ZZ). 



(6) 



The second term in 5i a tt originates from the measure D[i(Z,Z) = D(Z,Z) SDet(l — ZZ). Note that the functional 
integral for the two-point function (O(M)O(N)) is regularized by the boundary condition (||). 

What we have done is an exact rewriting of the original problem. It might superficially look as though we have 
made the problem more complicated by transforming from U(l) phase integrals to an integral over supermatrix fields 
Z and Z, but this is not so. The point is that while the U(l) phases Uq(1) fluctuate independently, the dominant 
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contributions to the integral over the supermatrix field come from slowly varying field configurations. The latter 
property is best seen by combining terms to rewrite the lattice action as follows: 



S latt [Z,Z] = In SDet H (1 - (1 - ZZ)- X Z{Z V - zj) . 



The right-hand side vanishes for constant fields (Zjj = Z) and is small for slowly varying ones. This will allow us to 
take a continuum limit and describe the low energy or long wave length physics of the Chalker-Coddington network 
model by a continuum field theory. In view of (Q) and , it should not come as a surprise that this field theory will 
turn out to be Efetov's nonlinear a model with unitary symmetry, augmented by Pruisken's topological term. 

The formula we have derived for {guN) illustrates how the transport coefficients of the network model can be 
expressed as correlations functions of a supersymmetric lattice field theory. In the remainder of the paper we will not 
discuss the specific correlator (<7mat) any further, but will concentrate on the general structure of the theory. 

V. GLOBAL SYMMETRY UNDER GL(2|2) 

To get ready for taking the continuum limit, we will now elucidate the symmetries of the lattice field theory. We 
are going to show that it is invariant under global transformations 

Z(l)^g-Z(l) = (AZ(l)+B)(CZ(l)+D)~\ 
Z(l)t-> g • Z(l) = (C + DZ(l))(A + BZ(l)y\ 

A B 



for 9 — i (j jj j S GL(2|2). The integration measure D (z(l), Z(l)^J has this invariance by definition. To see that 



the lattice action (Q) is invariant, we first transform the factor SDct ^1 — Z(l)Z(l)J for a single link I, temporarily 
dropping the argument I for notational simplicity: 



SDet (l-(g- Z)(g ■ Z)\ 

SDet (l - (AZ + B)(CZ + L>)" 1 (C + DZ)(A + BZy 1 



= SDet ( 1 - (Z + A^B^l + D~ 1 CZy 1 (Z + D^Cyi + A^BZy 1 



This expression is further processed by setting A B =: X and D C — : X and using the identity 

i-(z + x)(i + xzy x (z + x)(i + xzy 1 
= (t - xx) (i + zxy l {\ - zz){\ + xzy 1 , 

which follows from elementary algebra. We then get 

SDet (f - (g ■ Z){g ■ Z)\ = SDet ((1 - XX) (1 + ZX)-\l - ZZ){\ + XZy 1 ^ . 

The "dynamical" factor SDet^(l — ZjjZy 1 is transformed in an identical fashion. Because the transformation is 
global and U\ acts as the identity in superspace, we have 

(g ■ Z) v = U x {g- Z)uy = g ■ {U^U^) = g ■ Z a , 

i.e., the actions of g and U± commute. We thus obtain 

SDet n (1 - (g ■ Z) v (g ■ Z))^ = SDet H (f - (g ■ Z v ){g ■ Z)^ 
= SDet w ((1 - XX)(1 + Z v X)-\\ - ZuZ){\ + XZ'y 1 ^ 

By combining factors and using the multiplicativity of the superdeterminant, we arrive at 

S Utt [g -Z,g-Z] = 5i att [Z, Z] - lnSDet„(l + ZXy 1 + lnSDct w (f + ZuXy 1 . 







a) 



b) 



U 



R 



D 
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U 

R 



FIG. 3. The two types of elementary vertex of the Chalker-Coddington model. Type b) is obtained from type a) by a rotation 
through 7r/2. 



The last two terms on the right-hand side cancel each other by 

SDet«(l + ZuX) = SDet w (l + U-yZU^X) = SDet«(l + ZX), 



Ut (i.e. 



U^XU-t = 



X). This establishes the global 



since X is constant on links and therefore commutes with 
invariance £iatt[sf • Z,g ■ Z] = S\ M [Z, Z). 

The global GL(2|2) symmetry of £i a tt is broken by the boundary condition (Eh, of course. Such a symmetry breaking 
is needed for the mathematical consistency of the formulation, since the integral over the q = 0, or zero momentum, 
mode of the noncompact bosonic sector would otherwise be divergent. 

In Sec. IV we saw that the supersymmetric reformulation of the Chalker-Coddington model is defined over the same 
field space as Pruisken's model. What we have just learned is that these models not only share the field space, but 
also have the same internal symmetries. This is already a strong indication of their equivalence at the critical point 
p = 1/2 resp. a^j = 1/2. 



VI. REFORMULATION AS A VERTEX MODEL 



Further analysis is facilitated by the observation that the "Boltzmann weight" of the lattice field theory, 

W[Z, Z] = SDet w (l - Z Zu)- 1 SBct H (l - ZZ), 

factors to a large extent. This factorization is seen as follows. Consider the expression Z(l")U 1 (l" ,l')Z{l')Uy (/' ,/) 
and take I to be the link / = D, say, emanating in the downward direction from the node shown in Fig. 3a. The 
matrix element ?7 1 _1 (Z', D) is nonzero only if V is one of the two links, V = L or /' = R in the figure, that flow into the 
same node. In either case, the matrix element Ui(l", I') leads to the final state I" being one of only two links, namely 
I" = I = D (the one we started from), or else I" = U, the link emanating from the node in the upward direction. 
Similarly, if the initial state is I = U, the final state is I" = U or I" = D. These statements remain true (mutatis 
mutandi) if the node of Fig. 3a is replaced by any other node of the network. (Half of the nodes of the network are 
obviously equivalent to the one considered. The other half of them has links emanating in the horizontal direction 
and nodes flowing into it in the vertical direction, see Fig. 3b. These nodes become equivalent to the one of Fig. 3a 
after a rotation by ninety degrees.) 

Thus, ZU 1 ZUi 1 viewed as an operator in the space of states supported on links, connects only pairs of links, 
namely the two links that emanate from one and the same node. This means that the Boltzmann weight W is a 
product of factors, one for each node: 



W[Z,Z]= [] n(z{U n ),Z{D n ),Z{L n ),Z{R n )Y 



(7) 



nodes n 



where the assignment of link labels U,D,L,R is defined by Fig. 3. 
for a single node works out to be 



From the definition of U\ in Sec. Ill, the weight 
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1l(z(U),Z(D),Z(L),Z(R) 



SDet ( [1 - Z(U)Z(U)][1 - Z{D)Z{D)][l - Z{L)Z{L)][l - Z{R)Z(R)]) /2 (8) 



xSDet" 1 f 1 _ P Z ( U ) Z ( R ) ~ (! - V)Z{U)Z{L) Jp(T=p)Z{U)[Z{L) - Z(R)} \ 
V y/p(l-p)Z(D)[Z(L)-Z(R)} 1 - pZ(D)Z(L) - (1 - p)Z(D)Z(R) J ' 

A model of the kind (j?]) is called a vertex model in sta tistic al mechanics, and 7?. is called the R-matrix. The global 
symmetries of this R-matrix will be investigated in Sec. VIII. In Sec. [x| we will show how to pass from the integration 



over the fields Z, Z to a summation over superspin degrees of freedom. 

VII. CONTINUUM LIMIT 

For values of the parameter p close to 1/2, where the Chalker-Coddington model undergoes a quantum percolation 
transition, the correlation (or localization) length is very large, and we expect the lattice functional integral to be 
dominated by fields that vary slowly. Our goal in this section is to extract from (^) the continuum field theory 
governing these slowly varying modes. 

Recall the following facts: (i) Pruisken's model and the supersymmetric (SUSY) reformulation of the Chalker- 
Coddington model are defined over the very same complex field space, G/H = GL(2|2)/GL(1|1) x GL(1|1). (ii) Both 
the action functional of Pruisken's model, S con t, and the lattice action of the SUSY reformulated Chalker-Coddington 
model, Si a tt, are invariant under global GL(2|2) transformations, (iii) The metric tensor of the field space G/H, 

g = STr(l - ZZ) _1 dZ(l - ZZy x AZ = -STr(dQ) 2 /8, 

is invariant under GL(2|2) [i.e. under Q i— > gQg~ x with g £ GL(2|2)], and this is the only second-rank supersymmetric 
tensor on the homogeneous space G/H that has this invariance. When combined with a standard field theoretic power 
counting argument, these facts lead to the conclusion that the continuum limit of the lattice action S\ a tt must be a 
linear combination of the four terms (/i, v = x, y) in ([!]), which are induced from g. 

A further constraint comes from the spatial symmetries. At large scales (i.e. in the infrared limit), where details 
of the network structure are washed out, the Chalker-Coddington model acquires an invariance under rotations in 
space. From the four terms C^ v one can form only two linear combinations that are rotationally invariant. These are 
precisely the ones that figure in the expression (|l|) for Pruisken's Lagrangian (Pruisken et al. |1(J] arrived at them by 
essentially the same argument), namely C xx + C yy and C xy — C yx . We thus conclude that the continuum limit of the 
SUSY reformulated Chalker-Coddington model at criticality (p = 1/2) is Pruisken's model at the critical topological 

coupling 6 = it = 2ira x J , and some as yet unknown coupling ai°J : 



Siatt [smooth fields] = S c 



=• , xy =1/2 



The only thing that remains to be done then, is to work out the numerical value of the coupling constant o~ xx . In the 
following subsection we will argue that 

4°i = 1/2 (for p = 1/2). 

The small value of a x °J = 1/2 means that the theory is at strong coupling. (The weak-coupling limit, where Pruisken's 
derivation is valid, is a x °J 3> 1.) The superscript (0) alerts us to the fact that this is a "bare" value, obtained by a 
naive continuum limit, i.e. without taking into account possible renormalizations coming from modes with short wave 
lengths. 

In the following we pay no attention to the boundary conditions and assume the network to be infinitely extended. 

A. The coupling constant a^J 

To determine the value of a x °J for the SUSY reformulated Chalker-Coddington model, we first expand Si^tt to 
quadratic order in Z and Z , and then take the long wave length limit. The quadratic part of S'latt is conveniently 
obtained by expanding the expression for the R-matrix: 
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-2 In 11 (z{U), Z{D),Z(L), Z{R)\ 

= STz(z{U) [Z(U) - 2(1 - p)Z(L)] + [Z(R) - 2 P Z(U)] Z(R) 

+Z{D) [Z(D) - 2(1 - p)Z(R)] + [Z(L) - 2pZ(D)] Z{L)^j + 0(Z 2 Z 2 ). 

Now recall the meaning of the labels U, D, L, R defined by Fig. 3a, and let (n x ,n y ) £ 1? be the Cartesian coordinates 
of the node in that figure. Then, if a is the lattice constant of the network, 

Z(L) = Z(an x — a/2, an y ), Z(D) = Z[an xi an y — a/2), 
Z(Rl) = Z(an x + a/2, an y ), Z(U) — Z(an x , an y + a/2). 

Similar expressions hold for the other type of node, shown in Fig. 3b (rotate by it/ 2). Summing In TZ over all nodes 
of the lattice we get an expression for the quadratic part £2 of Si a tt, which is of the general form 

S 2 = J2 W - T ( Z ' STtZ{1)Z(1). 

IV 

The matrix T(l,l') represents a bistochastic process, i.e. it obeys the sum rules 

]Tt(u') = i = EW)- 

1 i' 

The first of these means that, if {Pt(l)} is a set of dynamical variables with dynamics Pt+i(l) = T(l,l')P t (l'), 
then the total probability is conserved: J^i PtiP) — const (independent of time t). It is not hard to guess the physical 
meaning of this: the dynamics at hand is nothing but the classical approximation to the full quantum dynamics of the 
Chalker-Coddington model. In other words, expanding the lattice action Siatt to second order in Z, Z is equivalent to 
neglecting all quantum interference effects and approximating the Chalker-Coddington model by its classical limit. 

The classical model has been looked at in a number of papers, most recently in Sec. I B of (|l), where it was shown 
that the dynamics at large wave lengths simplifies to classical diffusion, with the longitudinal conductivity of the 
network being a x °J = \/p(l — p)- (There is also a nonvanishing Hall conductivity a x °y = P-) By translating this result 
into the present context, we obtain for S2 a diffusive action, 

\ZpO--p) J STr (-Zd 2 x Z - Zd 2 y Z) d\ 

for large wave lengths. Matching this expression to the quadratic part of S'cont, Eq. (Q), gives ai°J = 1/2 at p = 1/2, 
as announced above. 



B. The coupling constant ai°J 

The critical points of the Chalker-Coddington model and Pruisken's model are at p = 1/2 and a x °J —1/2 (mod 1) 
respectively. Therefore, as was argued earlier, there cannot be any doubt that the continuum limit of the lattice 
action (^) at p — 1/2 is Pruisken's action with topological coupling a x °J — 1/2. Nevertheless, it is both reassuring 
and instructive to check this by direct calculation, which is what we do next. 



In Sec. IVII A| the value of a x °J was found by looking at the quadratic part of the action. The topological coupling 



a xy\ by its very nature, evades any such attempt at perturbative calculation. To extract it from (|6|) a different, 
nonperturbative scheme must be used. A direct approach would be to perform a gradient expansion around a 



topologically nontrivial background. For reasons that will be explained at the end of Sec. VIII. this is not easy to do. 
Here we will follow a different procedure, which is to evaluate Siatt on the lattice discretization Z^ of some smooth 
field configuration with topological charge m/0. 

Given the relation ImS , CO nt[Z < - m \ Z^] = 2itmo- x °y and the requirement S'cont = Si at t for smooth fields, the topo- 
logical coupling is determined by 



ImS lat t[Z (m) ,^ (m) ] = 27rm4°). 

a x °J from this equation, it is 

special form 



To calculate a x °J from this equation, it is easiest to consider fields that have topological charge m = 1 and are of the 
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Here the components .Zbb, Zbf, and Zfb, which are topologically trivial, have been set to zero and only Z-$y nas 
been retained. With this choice, the formula for a x ° y reduces to 

^ = ^ImSiatt^Wz* 1 )] = -i-ImInDet w (l + fUJU^ 1 ) 



p £ Imln7e FF (/([/„), f(D n )J(L n )J(R n )) 



27T 

nodes n 



where 7^ff is the R-matrix in the FF sector: 

^FF {f(U n )J(D n )J(L n ),f(R n )) 

not ( 1 + /(WW + /( i? )l/2 /(WW - /(«)]/2 \ 
" \ f(D)[f(L)-f(R)]/2 l + f(D)[f(L) + f(R)}/2j 

= 1 + \[f{U) + f(D)][f(L) + f(R)} + f(U)f(D)f(L)f(R). 
Now consider in the 2d plane with coordinates x and y the smooth held configuration 

f{x, y) = (^p- + l X -p zo) 1 (A ± e R), 

which interpolates between / = oo (corresponding to the south pole on the two-sphere Mp = S 2 ) at 

1 — i 1 + 2 
x + iy~ — - — A-Rcz -\ — A + Imz (9) 

and / = (the north pole on S 2 ) at infinity, and has topological charge m = 1, as can easily be checked. (Note that 
/ is not a onc-instanton configuration, unless A_ = A + .) The parameters A± set the size in x ± y of the region where 
the topological density Ltop for / is appreciably different from zero. 

For this choice of /, we now work out the R-matrix at the node in Fig. 3a with coordinates (n x ,n y ). Setting 
a± = a/(2A±), and C = [n x — n y )a/A_ + i(n x + n y )a/A + — z , we have 

f(U) = (C - a_ - ia + y\ /(£>) = (C + a_ + ia+J" 1 , 
/(L) = (C - a- - f(R) = (C + a- + 

The resulting value of the R-matrix is 

2|C| 2 + 1 



TZ FF = 1 



[C 2 - (a- + la+YWC - (a- + ia + f] 



In order for the field configuration to be smooth on the lattice, the parameters a± must be small. This allows us to 
Taylor expand with respect to one of these, say a+. Taking also the imaginary part of the logarithm, we obtain 

±mln"K F F = 4a+a_p(l + p) _1 Re (( 2 - o^)" 1 + 0{a\), 

where 1 + p = 7?-ff|q + =o- 

In the next step we sum the contributions from all nodes (n x + k,n y + k) with k e Z. The smallness of a + <C 1 
allows us to convert the expression for ^Imln7^FF into an integral. Decomposing £ into real and imaginary parts 
by C = £ + if], 

£ = 2a-(n x - n y ) - Rcz , V = 2a + (n x + n y ) - Imz , 

we arrive at 

Vlmln^FF / drj-^- Re^ + ^f-a^y 1 . 
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This integral is easily evaluated by closing the contour and applying the method of residues. One finds 



What is the geometric interpretation of this result? Let us agree that the word "vertex" here means a node 
taken together with its four links, severed at half the distance to neighboring nodes. The set of vertices with node 
coordinates (n x + k,n y + k) (k £ Z) sweep out a diagonal strip of the 2d network. The above result means that 
Imln7?.FF vanishes if the center @ of the topological excitation lies outside the diagonal strip swept out by the 
vertices (n x + k,n y + k) (with variable k £ Z), and equals — it when it lies inside. 

The above calculation applies to the type of node shown in Fig. 3a, which may be characterized by the condition 
n x + n y £ 2Z, say. Doing the calculation for the other type of node (n x + n y £ 2Z + 1 , or Fig. 3b) gives the same result 
except for a change of sign — 7r — > n. Now observe that the diagonal strips swept out by the vertices (n + k, —n + k) 
and (n + k, —n + 1 + k) (k,n £ Z) cover the plane completely without overlapping. Therefore, 

exp Tm5iatt[^ (1) ,^ (1) ] = exp ^ -TmlnftpF = e ±m = -1, 

all nodes 

independent of the location of the topological singularity. By a x ° y — (27r) _1 Im5iatt[^^ 1 \ Z^], this proves a x °J = 1/2 
(mod 1). 



C. Generalization to N channels 

D.K.K. Lee and Chalker p3] introduced a generalization of the network model that has two channels per link. The 

one-step time evolution operator of that model, U^ 2 \ is again a product of factors: = Jj[ Uq \ . The second 

factor is diagonal on links and associates with each link I a 2 x 2 matrix Uq{1) drawn at random from the unitary 
group in two-dimensional channel space, U(2). The first factor describes the deterministic scattering at the nodes. 
Depending on the choice made for this factor, the two-channel model applies to a spin degenerate Landau level or 
electrons in a random magnetic field (the so-called random flux problem). 

In this subsection we consider an iV-channel generalization of the one- and two-channel network models, where 
random U(iV) matrices are placed on the links. The one-step time evolution operator is written [/ W = Uq . The 
two factors describe the deterministic iV-channel scattering at the nodes, and the random U(iV) directed propagation 
along links, respectively. 

The case of N channels per link can be treated by a slight extension of our field theoretic formalism. Such an 
extension is possible since the basic formula @j was shown in JTHJ to generalize from U(l) to XJ(N), as follows: 



f dU exp U\ a U^^ +a + #. T ffV-r) 

JV(N) V 7 

D m (Z, Z) exp (i,\ a Z aT ^ T + ^ T Z TCT V>i ff ) , 



Efotov 



where U = Uq(1), and the link label / was suppressed for clarity. The right-hand side differs from that of (|J) only by 
the channel index i = 1, N attached to ip, tp, and the different form of the weight function in the superintegration 
measure: 



Dfx N (Z,Z) := D(Z,Z) SDet(l - ZZ) 



N 



Using the above generalization of (Q), we can reformulate the iV-channel model as a supersymmetric lattice field 
theory with action 



S^[Z, Z] = lnSDet W(5C iv 1 - Uf'ZUf n Z) - JVlnSDet«(l - ZZ). 



The derivation exactly parallels that of Sec. IV. The first superdeterminant runs over the tensor product of superspace 
with link space (Tt) and channel space (C^ ). 

If we choose the deterministic scattering u[ N ^ to be of the special form — Ui~^ ® ljy, i-e. u[ N ^ acts as the 

identity in channel space, is simply a multiple of the action of the one-channel model: 
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S^ t [Z,Z]=N,S^[Z,Z]. 



Therefore, from Sees. VII A and VII B the coupling constants of the corresponding continuum field theory are 



<&> = = N/2. 

The choice u[ = U^ N ® 1 at for N = 2 is appropriate |22] for the spin degenerate Landau level at the center of the 
band, and for the random flux problem at the symmetric point. Because Pruisken's model at a^J =2/2 (mod 1) = 
is known to be a massive theory, our result for a^y confirms the proposition of j2^] that neither of these systems is 
critical. Rather, they are in a (Haldane type) massive phase [Q, corresponding to localization of all states. 

VIII. SYMMETRIES OF THE R-MATRIX 

The basic building block of the SUSY reformulated Chalker-Coddington model is its R-matrix. For the case N = 1, 
which is to be analysed here in more detail, this building block was given in (||). To gain a deeper understanding of 
the model and possible variations thereof, we are now going to investigate the symmetries of that R-matrix. Most of 
the effort will be expended on rewriting the expression (JsJ) in a different form, so as to make those symmetries more 
evident. 

We start by undoing the tp, ip integration to write the R-matrix as 

11 = II SDct i 1 - Z(0)Z(0)^j 1/2 SDet (l - Z{I)Z{I)\ 1/2 
i,o 

x J D(V>,$exp(- i> +a (0)iP +<T (0) +$ +a {0)Ui (O, J)Z OT (/)^_ r (J) 

-^_ T (7)^-r(I) + 4>- T {I)U^{I,0)Z Ta {0)^ +a {0) 

Here we have introduced the labels I £ {il, i2] := {L, R} and O £ {ol, o2} := {U, D}. The notation is motivated by 
observing that, according to the direction of motion indicated by the arrows in Fig. 3, the links L and R are incoming 
states for the scattering at the nodes, while the links U and D are outgoing states. 

To proceed, we need to recall briefly various mathematical structures that were developed in detail in the appendices 
of jl5| and p3[ . First of all, we introduce Fock operators c and c, which are quantum counterparts of the classical 
variables and ip. Let b±,b± and f±,f± be canonical boson and fermion creation and annihilation operators, and 
set 

c+f = /+, c +B = b+, c_f=/I, c_b=&1, 

C+F = /+, C +B = 6+ , C_ F = /-, C_ B = -&— 

The operators c and c are canonical pairs, with graded (or super) commutation relations 

[cx,Cy] := C X C Y - (-l) |X||r| CFCJC = S X Y, 

where \X \ = if X = ±B and \X\ = 1 if X = ±F. They act in a Bose- Fermi Fock space with vacuum 

c +B |0) = c +F |0) = c_ B |0) = c_ F |0) = 0. (10) 

The graded commutation relations are invariant under canonical transformations, cx l— * ^g^x^g 1 ~ ^y9yx and 
c x h-> TgCxTg 1 = (g- 1 )xYC Y , where 



C £> 



G GL(2|2) i ^ T s := exp (c x (lng)xYC Y ) 



defines a representation of GL(2|2) on Fock space. 

Consider now the subspace, V , selected by the condition excx — (summation convention!) or, equivalently, by 

b{b + + flf + = blb_ + /!/_. 
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The resulting constraint on the Bose- Fermi occupation numbers is rib+ + n s+ = n b- + n f-- By using this equation 
to eliminate rib-, say, we can characterize the states of V by a triplet of integers (ri6+, ^/+, n f-)i where nf + and 
take values from the set {0, 1}, and rib + = 0, 1, 2, oo. 

Let Py be the operator that projects Fock space onto the subspace V . As was shown in (TH), this projector has a 
resolution 

Pv = J D{Z,Z) \Z)(Z\ 

by generalized coherent states, 

\Z) = exp(c +(T Z . T c„ T ) |0) SDct(l - ZZ) 1/2 , 
(Z\ = SDet(l - ZZ) 1/2 (0| exp f-c_ T Z T(T c +(T 



By simple manipulations exploiting the standard properties of coherent states, one can verify the following equality: 

Y[ SDet (l - Z(I)Z(I)) V2 x exp (^(0)^(0, /)Z ffT (/)^_ T (/)) 

= (0|exp(^ +CT (o//)J7i(o/x,ii/)c +CT (i/)-c_ T (j/)V'-r(w)) \Z(il) <g> Z{i2)) 
= (0| exp (^ +(T (o/i)c +CT (/i) - c_ T (z/)V>- T (w/)) 

x exp fc +ff (M)(ln^i)(M, f)c+aO)) |Z(il) <8> Z(i2)>. 



Here /i, f € {1, 2}, and C/i is defined by identifying initial and final channels, i.e. Ui(fj,, v) — Ui(ofi, iv). Similarly, 

1/2 

J] SDet (l - Z(0)Z(0)) x exp (^_ T (7)C/f 1 (/, O) Z TtT {0)^ +(7 {0)) 
o 

= (Z(ol) (g) Z(o2)| exp (i/L T (ii/)J7f 1 (w, o/i)c_ r (^i) + c +a (ti)i> +IT {ofi)) |0) 
= (Z(ol) ® Z(o2)| exp (c_ T (z;) (In #i)(^/x)c_ t (m)) 

x exp (c +<T (fj,)ip +<7 (ofi) + i/)_ T (ii/)c_ T (i/)) |0). 

The variables V'jV' have now served their purpose and we integrate them out, by using the closure relation for Bose- 
Fermi coherent states: 



id = J D(i/j, ip) exp (-ip +a (v)ip+a(ti) - ijj- T {v)i>- T (v)) 

x exp (c +rT {pb)ijj +rT {pb) + i)- r {v)c- T {v)) |0) (0| exp (ip +a (n)c +<T (n) - 5- T (u)ip- T (p)) . 
All this results in the following formula for the R-matrix: 

K (z(U),Z(D), Z(L), Z(i?)) = (Z(U) <g> Z{D)\n\Z(L) ® Z(R)}, 
where the operator 1Z is expressed by 

K = exp(c x (ti)(lnUi)([i,p)cx(v)) ■ (11) 
The advantage of this reformulation is that the invariance of 1Z under GL(2|2) transformations 

cjc(m) ^ TgCxi^Tg 1 = c Y (fi)g Y x, 



is obvious whereas previously, in Sec. [Vfl we had to work quite hard to establish the global GL(2|2) invariance. 
Although 1Z commutes with T g , its matrix elements are not invariant: 



K Z 



!(U),Z(D), Z(L), Z{R)) ?1l(g- Z(U), g ■ Z(D),g ■ Z(L),g ■ Z(R)) 



1G 



It is worth spending a little effort to explain how that comes about. As will be seen, the reason is that, since the 
Fock space vacuum is not a scalar with respect to H — GL(1|1) x GL(1|1), the coherent states do not transform as 
functions on G/H, but rather as sections of an associated line bundle, see (l| and references therein. In other words, 

T g \Z) + \g-Z). 

The correct transformation law is derived as follows. We write the coherent states as \Z) = T s ^ z jh|0), where 



s(Z,Z) 



{l-zz)- 1 ' 2 zix-zzy 1 ' 2 
z(i-zz)- 1 / 2 {i-zzy 1 / 2 

1 Z\ ( {I- ZZ) +1 I 2 \/l 

1 ) \ o {i-zzy 1 ' 2 ) \z 1 



We then define an 77-valued function h(g; Z, Z) by 

9 s{Z, Z) = s(g ■ Z,g ■ Z) h(g; Z, Z). 

The explicit form of h(g;Z 7 Z) can be found in Appendix B of [jl3). From ( |l0| ) and the definition T g 
exp (cx(ln g)xYCY), one easily sees that the vacuum carries a one-dimensional representation fj, of H: 

T h \0) = |0) fj,(h) (for he H), n (diag(A, D)) = SDetZ)- 1 . 

Therefore, the coherent states transform as 

T 9 \Z) = T g T s{zrz) \Q) = T s{g . Zg . z) T hig . zz) \0) = \g ■ Z) /x (h(g;Z,Z)) . 

As a result, the R-matrix obeys the following transformation law: 

TZ (g ■ Z(ol), g ■ Z(o2), g ■ Z(il), g ■ Z(i2)) = K (z(ol), Z(o2), Z(il), Z{i2)) 

\{ ti(h(g;Z(o\),Z(oA))) J[ (i (h(g; Z{iv), Z(iv))^J _1 . 



x 

A=1.2 t/=l,2 



Because each link is incoming with respect to one node, and outgoing with respect to another, the multipliers 
(i (h(g; Z, Z)j cancel when all R-matrices are multiplied together, so that the global GL(2|2) invariance is recovered. 

The transformation law for the R-matrix explains why a direct gradient expansion to extract axy from Siatt is 
difficult. Such an expansion locally produces terms such as 

STr(l - zzy x zd x z(\ - zzy x zd y z, 

for example, which are singular at Zpp — oo (the south pole on Mp = S 2 ). When all of these terms are correctly 
summed over the entire network, they cancel, as is guaranteed by the global GL(2|2) invariance, which permits to 
rotate the south pole into any other point on the two-sphere S 2 . However, the cancellation really does take place only 
after summation of terms. By the multiplier-corrected transformation law of the R-matrix, singular terms remain 
locally, making the extract ion of the topological coupling difficult. This, then, is the reason why a gradient expansion 



was not attempted in Sec. VII B 



IX. ANISOTROPIC LIMIT 



We have presented an analytical method for dealing with the Chalker-Coddington model in its original isotropic 
formulation, by mapping it on a lattice equivalent of Pruisken's nonlinear a model. In this section we will review 
another way of arriving at Pruisken's model, a replica version of which was first published by D.H. Lee. Following Q] 
we now take for our starting point the anisotropic limit of the Chalker-Coddington model, and replace the unitary 
operator U = U\Uq by the Hamiltonian H for an array of chiral modes n = 1,2, ... with velocity v and alternating 
direction of propagation: 
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H = jdx^ *Ux) [5 nn ,(-l) n ivd x + V nn ,{x)] 



(x). 



The functions V nn i(x) = V n i n (x) are uncorrelated Gaussian random variables with zero mean and variance 

(Vnn'(x)V nn >(x')) = 2 (u S nn > + Ui6 rhn > + 1 + Ul5„, n '-l) K x ~ x ')- 

The symbol § means that we are using periodic boundary conditions in x. 

To prepare the treatment of the general case, we shall first consider the case of a single chiral mode n. The 
supersymmetric generating functional for the retarded and advanced Green's functions of H is set up in the usual 
way, see Sec. [V. Ensemble averaging over the random potential V(x) — V nn (x) leads to the functional integral 

Z = I V(tp,tp) exp f dx ^jj(Avd x - e)ijj - u (V>AV>) 2 ) , 



where s is a positive infinitesimal. As before, tpx (x) is a super "spinor" field with four components X = +B (retarded 
Boson), X = +F (retarded Fermion), X = — B (advanced Boson), and X — — F (advanced Fermion). The notation 
means ipAijj — 4>+ a ip+ a — ip-rip^r- If the energy in the retarded and advanced sectors is different, lu = E + — E_ ^= 0, 
we need to add a term iujipip to the Lagrangian. To probe this field theory more generally, we may couple it to an 
external nonabelian gauge field A{x) S Lie(GL(2|2)) and consider 

Z^irM := J V(rP,ij) exp I dx {^Av{d x + A)i> - u(^A^) 2 ) . 

We have temporarily set u = uq for notational simplicity. The special coupling to frequency is retrieved by putting 
A = (iuj — e)A/v, independent of x. Note that the generating functional -Z^ irac [yl] is invariant under local gauge 
transformations, 

Zmr,M = 2^ c [ h A], h A = hAh- 1 + hd x h-\ 

where h(x) S GL(2|2) acts on the spinor field by ijj i— > hip and i— > ipAh -1 . 

It turns out that one can write down another 0+1 dimensional field theory that has the very same local gauge 
invariance. The field of this theory is the supermatrix Q = gAg -1 , which was defined in Sec. || and transforms as 
Q i—* hQh^ 1 . The generating functional is 



Z WZ [A] := J VQ exp I dx STr ^-cT 1 ^ + A)g. 



Clearly, this satisfies Zwz[A] = Zwz[ h A\. The linear derivative term in the action is of the Wess-Zumino type, i.e. 
cannot be expressed in a globally nonsingular way in terms of Q, and is often called a Berry phase. The theory 
is well-defined because the ambiguity under right translations g ^ gh^ [h^Ah^ 1 — A or, equivalently, /ir(x) 6 
GL(1|1) x GL(1|1)] gives rise to a factor 



exp (J) dxSTr—h^dxhn = exp j) dxd x STr— In /ir = exp 2 



nim = 1, 

which is unobservable in the functional integral. (For more details see Sec. 3.4 of [fl3f .) 

The local gauge invariance shared by -2^.^, [^4.] and Zwz[A] suggests the existence of some relation between these 
theories. In fact, the following statement is true: 

lim Z£ iiac [A] = Z WZ [A}. (12) 

u — >oo 

This identity can be viewed as a + 1 dimensional analog of nonabelian bosonization in 1 + 1 dimensions and, since 
A couples to vip'tpA and gAg~ 1 /2 = Q/2 in the respective cases, amounts to the "bosonization rule" 

v^A ™ Q/2. 

We now briefly sketch the proof of the nonabelian bosonization formula ( |l2| ) . By the local gauge invariance of both 
theories, it is sufficient to prove the equality for an ir-independent gauge field A. Moreover, A may be taken to be 
a diagonal matrix. In this special case, it is easy to apply a method that was described at length in p3] and works 
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as follows. As a first step, one identifies -2-p irac L4] as the coherent state path integral of a supersymmetric Hubbard 
Hamiltonian for bosons and fermions. Then, one takes advantage of the limit u — > oo, which enforces a Hubbard 
constraint reducing the (low energy) degrees of freedom to that of a single superspin. And finally, one sets up the 
coherent state path integral for the superspin Hamiltonian. The latter path integral turns out to be Zwz[^4], which 
completes the proof. 

This proof, although straightforward, has the disadvantage of being somewhat indirect. A more direct procedure is 
to decouple the interaction term (ipAip) 2 by introducing a Hubbard-Stratonovitch field Q coupling to iptpA and then 
to integrate out ip and ip. The effective action for Q is 



S[Q] = jdx STr (-^Q 2 + Hdx + A + uQ/v 2 )^) . 



The next step is to simplify the Q field functional integral by means of the saddle-point approximation, as a result 
of which Q gets restricted to the nonlinear space Q = gAg -1 . This step, while only approximate in general, here 
becomes exact in the limit u —> oo. [What makes this possible is the stationarity of the average density of states of 
vid x + V(x).] By expanding In [Ag^ 1 (d x + A)g + u/v 2 ) to linear order in g~ 1 (d x + A)g, one obtains the action of the 
Wess-Zumino functional Zwz[^4]- Higher orders are suppressed by powers of v 2 /(L x u), with L x the system size. 

Let us finally return to the case of many counterpropagating chiral modes that are coupled by hopping matrix 
elements between neighboring modes, with variance u-y. The Gaussian random hopping gives rise to an additional 
term in the Lagrangian, 

L L + 2u% ^ (V'nAVvj+i) (^ n+1 Aip n ) 

n 

= L + 2u\ 22 STr (V'nV'nA) (^n+iVVi+iA) . 
n 

By the bosonization rule vip'ipA — ► Q/2 for uq — ► oo, the additional term turns into (ui/2v 2 ) ^2 n STr(Q n Q n+ i) . The 
condition of validity of this step is uq ^S> u\. As a result we obtain the Q field action 



S[Q] = idxY, STr ((-l)n±g-ld x9n + ^QnQn, 



By a standard calculation [g(j,^3| this is the action of the coherent state path integral for a quantum superspin 
Hamiltonian, 

V n XY 



where S XY — excy, and the graded commutation relations of the Fock operators c, c were given in Sec. VIIL To 
reproduce the alternating sign of the Wess-Zumino term, we must alternate the definition of the Fock vacuum. On 
even sites (n € 2Z) the relations (|h]) apply, whereas on odd sites (n £22+1) we have 

c +B |0) = c +F |0) = c_ B |0) = c_ F |0) = (14) 

instead. The alternating vacuum plays the same role as the Neel state for ordinary spin systems and means that the 
superspin chain is "antiferromagnetic" in character. 

The Hamiltonian H spin with translational invariant coupling J = 2u\/v 2 > was shown in [T^ ] to represent the 
low energy limit of the quantum Hamiltonian of Pruisken's nonlinear a model at criticality. Thus, we finally conclude 
that the anisotropic Chalker-Coddington model with homogeneous (on average) inter-channel hopping is in the same 

universality class as Pruisken's model at o™ = 1/2. (Clearly, this line of reasoning is much less direct than that 



presented for the isotropic model in Sees. IV VII.) 

To conclude this section let me mention that the nonabelian bosonization formula ( |l2| ) extends to N channels: 



J VQ exp j> dx NSTi^g- 1 ^ + A)g 



= Jirn^ J V(ip, ip) exp j> dx (ip n A(d x + A)ip n - u{ip n A%p n ,){il) n , Aip n )) . 

The strategy of the proof is the same as for N = 1. Using this formula we can easily reproduce the A^-channel 
result a x °J — N/2 of Sec. VII C. The iV-channel bosonization formula offers also a quick way of analysing the 



zero-dimensional limit of the 2d chiral metal 
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X. MODIFIED NETWORK MODEL 



The primary goal of all field theoretic analysis of the plateau-to-plateau transition in integer quantum Hall systems 
must be to identify the fixed point theory that describes this transition and uncover the conformal structure it is 
expected to have. Recent attempts p3 in this direction started from the observation that Pruisken's model or, 
rather, the closely related Dirac theory with random mass, random scalar potential, and random gauge field, has a 
global GL(2|2) symmetry for the case of one retarded and one advanced Green's function. (A finite imaginary part 



of the energy argument of the Green's functions reduces this symmetry to GL(1|1) x GL(1|1). However, Sees. [II and 



IV 



show how such a symmetry breaking can be avoided by calculating a conductance between interior contacts of the 
network model.) This symmetry was then assumed to be promoted in the infrared to a Kac-Moody symmetry, which 
severely restricts the number of possible candidates for the fixed point theory. Unfortunately, these attempts have 
not been successful so far. What is needed as additional input to such considerations, which are purely algebraic, is a 
firm understanding of the Hilbert space structure, or the representations involved. It is one of the aims of the present 
work to contribute to such an understanding. 



Following up on unpublished work by N. Read, it was argued in 13 that the quantum Hamiltonian of the critical 
theory should be a superspin Hamiltonian of the type ( |l3| ) acting on a space of states built from alternating GL(2|2) 
modules 

... ® V V* ® V ® V* <g) 

where V and V* are generated by the action g i— > T g = exp {cx (In g)xYCy) of GL(2|2) on the vacua ( |To| ) and jl^ ) 
respectively. (As follows from Sec. 4.4 of jli|, the elements of V and V* have an interpretation as holomorphic and 
antiholomorphic sections of a line bundle associated to G — > G/H by the one-dimensional representation fj, of H. This 
permits the construction of a nondegencratc pairing between V and V* , so that these spaces can be viewed as being 
dual to each other, as suggested by our notation.) It is then natural to ask whether one might be able to construct an 
integrable superspin Hamiltonian, offering the possibility of an analytical and exact computation of critical properties. 



IX 



This question will now be addressed in the light of the results of Sees. VIII and 

Recall the "functional" vertex model presentation (0J|) of the Chalker-Coddington model, 

Z = Jv(Z,Z) II R(z(u n ),Z(d n ),Z(l n ),Z(r n j), 

nodes n 

R (z(l), Z(2), Z(3), Z(4)) = (Z(l) ® Z(2)\R\Z(3) ® Z(4)), 

R = exp (cx(jJ>)QnUi)(ji,, v)c x {v)^j ■ (15) 

(Note the change in the notation for links from capital to small letters.) It is instructive to pass from the integration 
over fields Z, Z to a summation over discrete degrees of freedom, as follows. Every link I emanates from one node, and 
ends at o ne node. Therefore, each Z(l) occurs once in the "bra" and once in the "ket" of some R-matrix. Recall from 



Sec. VIII that Py denotes the projector from Bose-Fermi Fock space onto the GL(2|2) module V, where the elements 
of V, referred to as superspin states, are labeled by three quantum numbers n/± € {0, 1} and nb+ = 0, 1, 2, ...oo. On 
using the closure relation 



f [] D ( z ( l ),Z(l)) <E> l \Z(l))(Z(l)\=<E> l P v (l), 



links / 

the lattice functional integral over Z, Z turns into a partition sum over superspin configurations {a(l)} [with a(l) £ 
V(l)] of a vertex model defined over the tensor product space ®iV(l): 

Z = ^2 (~1) Nf Yl R a(u n ) a (d n )Ml n )»(r n ) > Raf3,j8 = (C< ® f3\Rh ® S) , 

{a(l)} nodes n 

where (—1)^ is a sign factor due to supersymmetry (the partition sum is a supertrace), and R : V(l) ® V(r) — > 
V(u) ® V(d) is still given by jl^). For p — (left turns only), R can be seen to be the identity map, while for p = 1 
(right turns only) we have R = V, where 

V\a®L3) = (- 1) H 1/31 1/3 <8> a) 
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a) 



b) 




FIG. 4. a) Illustration of the row-to-row transfer matrix of the supersymmetric vertex model associated with the one-channel 
Chalker-Coddington model, b) The transfer matrix of the anisotropic model at p = translates the system by one lattice unit. 



is the graded permutation operator. 

Let me mention in passing that the presentation as a superspin partition sum can also be obtained from the network 
model directly p5f , without passing through the intermediate stage of a Z field formulation. 

How is the model built from the vertex R a /3,-y8 related to the superspin Hamiltonian (p"3|)? Given R a /3^s we can set 
up the row-to-row transfer matrix, T, illustrated in Fig. 4a, by summing over the superspin degrees of freedom that 
are situated on the horizontal links. (The meaning of the arrows in the present context will become clear below.) The 
corresponding Hamiltonian, defined as the logarithm of the transfer matrix, is nonlocal in general. A local Hamiltonian 
results on making the following modification of the isotropic network model. On nodes with coordinates (n x , n y ) G Z 2 



such that 



€ 2Z, the R-matrix is taken as it stands; but on the other half of the nodes (n x + n y € 2Z + 1), we 



replace the left-right asymmetry parameter p by its complement 1 — p. (At the level of the random network model, 



this is precisely what one does to arrive at the anisotropic limit of Sec. IX.) The row-to-row transfer matrix of the 
resulting anisotropic vertex model has the property, for p = 0, of translating the system by one lattice unit (Fig. 4b). 
By differentiating the logarithm of this transfer matrix at p = 0, one gets a superspin Hamiltonian that couples only 
nearest neighbors and is precisely the Hamiltonian H sp i n of (13). 



f/spin is not expected to be exactly solvable. However, one may ask whether it could be made so by slightly 
changing some parameters while keeping the general structure and symmetries the same. To get some hint, we turn 
to the well-developed theory of integrable systems p7fl . There, the integrability of a Id quantum Hamiltonian H 
is traced back to the existence of a transfer matrix T(u) that depends on a "spectral" parameter u in such a way 
that T(u)T(v) = T(v)T(u) for all u, v, and H is the logarithmic derivative of T(u) at some special point u = u$. A 
sufficient condition for T{u) to form a commuting family is known to be the quantum Yang-Baxter equation for the 
R-matrix. 

To my knowledge, most (if not all) of the integrable models discussed in the literature have one characteristic feature 
in common: their commuting family of transfer matrices possesses a "classical" limit, where T becomes the identity. 
In contrast, for the row-to-row transfer matrix associated to the Chalker-Coddington model (Fig. 4a), no such limit 
exists. The reason is simply this. The R-matrix R :V ®V — > V ®V hi ^) relates incoming channels on horizontal 
(or vertical) links to outgoing channels on vertical (or horizontal) links. On the other hand, the transfer matrix T 
propagates the degrees of freedom from one row to the next. Therefore, to construct the row-to-row transfer matrix 
from the R-matrix, we must reinterpret some initial states as final states, and vice versa. This is done by noting that 
the space of linear maps V — - > V is isomorphic to the tensor product V <g)V* . In this way, one sees that in Fig. 4a 
vertical links with an arrow pointing up carry the space V, whereas vertical links with an arrow pointing down carry 
the dual space V*. (This, then, is the meaning of the arrows in that figure.) Hence, the row-to-row transfer matrix 
of the SUSY reformulated Chalker-Coddington model is a map 



T : 



V ®V* 



V < 



V* <8> V <g> V* 



which connects inequivalent spaces. Therefore, T cannot ever be the identity. This means that the Chalker-Coddington 
model lies outside the category of vertex models for which the well-developed Yang-Baxter machinery applies. Thus 
it seems that there exists no known systematic way of deforming the Chalker-Coddington model to integrability and 
obtain an analytical solution. 

The above discussion, though being a falsification, also suggests a remedy. Given that the standard formalism of 
the theory of integrable systems requires the row-to-row transfer matrix to be a map 



T : 



V ® V* O V < 



V ® V* ® V < 



(16) 



we can turn things around and modify the network model accordingly. Consider the network shown in Fig. 5a. 
As before, the electron follows the direction of motion dictated by the arrows. It picks up a random U(l) phase 
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a) 



b) 



FIG. 5. a) One-channel Chalker-Coddington model, modified so as to make the direction of motion invariant along every 
horizontal and vertical line, b) The model of a) modified further, by doubling the number of channels on the horizontal links. 

while propagating along the links, and is scattered deterministically at the nodes, just as for the original one-channel 
Chalker-Coddington model. The crucial difference from Fig. 2 is that now the direction of motion does not alternate 
but is invariant along vertical and horizontal lines. An electron incident on a node either passes straight through it, 
with probability p, say, or else is scattered to the right or left, as the case may be, with probability 1 — p. (Note 
that such a modification of the scattering dynamics has no justification from a microscopic picture of guiding center 
drift along equipotentials. However, since our aim is only to describe the critical behavior, we should have a certain 
amount of freedom in the choice of model.) 

By construction, the row-to-row transfer matrix of the supersymmetric vertex model associated with the network 
model of Fig. 5a is a map of the desired type ([jl]). One may now hope to be able to deform and extend this 
transfer matrix to a one-parameter family of commuting transfer matrices. In somewhat more detail, this hope 
is based on the following facts. Consider the Lie superalgebra Q of polynomial maps u G C — > gl(2, 2) (with u 
being the spectral parameter). If Eab are the canonical generators of gl(2,2), the classical r-matrix r(u,v) = 
(u — v)^ 1 J2ab Eab{—^ B \Eba gives rise to a co(super)commutator Q — > Q®Q in the usual way fl6|| , thereby turning 
Q into a Lie bisuperalgebra. Quantization of Q leads to a Yangian superalgebra y :— Y (gl(2, 2)), which is a ^-graded 
Hopf algebra deformation of the universal enveloping algebra of gl(2, 2). Let A : y — > y (g) y be the comultiplication of 
y, and A op = Vo A its opposite. According to general principles p8| ], there exists a formal object called the universal 
R-matrix 1Z G y <£> y, which is determined by the intertwining relation A op (a) = TZA(a)TZ~ 1 and has the expansion 

InK(u) = u- 1 Y2(- 1 ) IBIe ab <8> E BA + 0(w~ 3 ). 

AB 

An irreducible matrix representation p of y yields an i?-matrix R p (u) = (p <S> p)(lZ(u)^ that is a rational function of 
u and solves the quantum Yang-Baxter equation. 

When applying this formalism to our problem, we should beware of potential problems due to the infinite- 
dimensionality of the spaces V and V*. Nevertheless, it does not seem unreasonable to expect the existence of 
R-matrices 

R(u) = (p v ® p v )(n(u)) : V ®V ® V, 
R(u) = (p v ® p v *)(K{u)) : V®V* -» V®V*. 

I have been able to verify that this expectation is fulfilled in the case of R(u), by explicit construction. The existence 
of R(u) remains an open question at the present time. (Help from experts on quantum groups would be very much 
appreciated.) The latter case is more complicated to treat because the decomposition of the tensor product V <X> V* 
into gl(2, 2)-irreducible subspaces involves a continuous series of representations of gl(2, 2) (see Sec. 5.2 of p~3[] ) , while 
in V <8> V only a discrete series appears. 
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Suppose now that both R(u) and R{u) exist, at least in some domain of the spectral parameter u. Then, since the 
quantum Yang-Baxter equation is automatically fulfilled, we can build a one-parameter family of commuting transfer 
matrices. The next question is: given the modifications we have made, is the physics of such a model still that of the 
plateau transition? Surely, obtaining an exact solution remains a far goal. It is therefore helpful that the mapping 
onto Pruisken's nonlinear a model provides us with a quick way to get oriented in the enlarged landscape of modified 
network models. Consider the model of Fig. 5a with parameter p — 1/2. Taking the continuum limit and computing 

the topological coupling in the same way as in Sec. VII B , one finds <Jxy = (mod 1). This means that the model 
does not lie in the quantum Hall universality class, but is in a massive phase with a finite (albeit large) localization 
length and exponentially decaying correlations. 

This result could have been anticipated from the following heuristic argument. Imagine separating the network into 
two independent subsystems A and B, one consisting of the vertical lines and the other of the horizontal ones. Then 
couple the modes within each subsystem by weak tunneling amplitudes. What you get in this way are two copies of 
D.H. Lee's anisotropic limit of the Chalkcr-Coddington model. By the reasoning reviewed in Sec. IX, each of these is 
critical, with the a model topological coupling constants being A = a^j B = 1/2. Now join the two subsystems 
to form the network model of Fig. 5a. From the meaning of Ltop as a topological density, it is reasonable to expect 
that &xy — a^y A + afy B — 2x 1/2 = (mod 1) for the coupled system, if the fusion is done in such a way that a^x 
for the coupled system is spatially homogeneous. Thus the network model of Fig. 5a will be noncritical at p = 1/2. 
(Its correlation functions at large scales should be similar to those of the two-channel Chalker-Coddington model at 
the symmetric point.) Since noncriticality is a generic property, this will remain so in a neighborhood of the point 
P=l/2- 

In view of this heuristic argument, we expect that criticality can be restored by superimposing on the noncritical 
network of Fig. 5a yet another copy of the anisotropic network model, thereby producing the network of Fig. 5b. There, 
the horizontal lines carry two, rather than one, channels per link. The propagation on vertical links is governed by 
random U(l) phases, as before, but the horizontal links now carry random U(2) matrices (just like the two-channel 
Chalker-Coddington model). One of the horizontal channels passes straight through the nodes, the other one is 
subject to the rules specified for the model of Fig. 5a. The topological coupling of the a model is then found to have 
the critical value a x ° v =3x1/2 (mod 1) = 1/2. Thus the last, doubly modified network model is critical and, on 
symmetry grounds, lies in the quantum Hall universality class. From what was said above, it is also a suitable starting 
point for attempting to deform to an integrable model. 



XI. SUMMARY 



Several messages result from the present paper. First of all, a close relation between two standard models of the 
integer quantum Hall transition, namely the Chalker-Coddington model at its symmetric point p — 1/2, and the 

supersymmetric formulation of Pruisken's nonlinear a model at 9 = — tt (mod 2ir), was established. Let us put 

this result in the proper context. To be sure, it has been clear for a number of years now that some sort of relation 
between these models ought to exist. We know that both are critical and belong to the same universality class, so 
they cannot but describe the same physics at long wave lengths. However, prior to our work, the understanding of 
the precise connection between Chalker-Coddington and Pruisken was rather indirect, relying on a double use of the 
anisotropic (or Hamiltonian) limit. The connection went as follows. At one end, by taking the Chalker-Coddington 
model and going to its anisotropic limit, D.H. Lee arrived at a network model consisting of an array of chiral modes 
with alternating direction of motion. At the other end, the anisotropic limit of Pruisken's model was investigated, by 
an elaboration of the work of Shankar and Read [^9| on the 0(3) nonlinear a model. It was argued in Jl3| that the 
Hamiltonian limit of Pruisken's model at 8 = tt and strong coupling (small a x °x) is an antiferromagnetic superspin 
chain. Now, the array of chiral modes and the superspin chain are easy to relate by conventional techniques. At 
the level of the replica trick this was done in the correct supersymmetric extension follows from (23). In more 
detail, the functional integral representation of the array of chiral modes maps on a nonlinear a model with an 
alternating sum of Wess-Zumino terms. (As we have seen, this mapping is based on a + 1 dimensional analog of 
nonabelian bosonization in 1 + 1 dimensions.) The latter, in turn, coincides with the coherent-state path integral of 
the antiferromagnetic superspin chain. 

One disadvantage of the above way of relating the models is that one does not have good control over the numerical 
value of the coupling constant a X x ■ In the present work, this uncertainty was resolved by dealing directly with the 
isotropic models, using a novel scheme devised in [ [l5|| . What we have shown is this. Starting from the Chalker- 
Coddington model at p — 1/2 and doing no more than an exact transformation followed by a continuum limit, we 
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arrive at Pruisken's model with 8 — n and a xx = 1/2. Or, in different words, the former model is equivalent to a 
specific lattice discretization of the latter. 

Previously, Pruisken's model was thought to be associated primarily with the white noise limit l c <C Ib, which is 
where Pruisken's derivation starts from. Recall, though, that in the course of deriving his model, Pruisken made a 
saddle-point approximation to eliminate the "massive" modes. Naive use of this approximation scheme is justified 
only in the limit of large o~ x °J (high Landau level). In contrast, the present work makes no such approximation. The 
only assumption we needed was the dominance of slowly varying fields in the functional integral, allowing us to pass 
from the lattice to the continuum. Thus, contrary to what might have been expected, Pruisken's model at small o~ x °J 
is actually associated more closely with the high- field limit, l c ^> Ib, as it is this limit that provides the microscopic 
justification of the network model. Note, however, that the ratio of microscopic length scales I c /Ib is expected to be 
an irrelevant parameter at a critical point with infinite correlation length. Thus, our result is not in conflict with 
Pruisken's derivation of the a model model as a critical theory. 

In my opinion, neither the Chalker-Coddington model nor Pruisken's model hold much promise for an exact ana- 
lytical solution in the near future. If so, the mapping of one model on the other is not yet a big step forward. The 
good news is that there are several useful spinoffs. D.K.K. Lee and Chalker suggested to model the random flux 
problem (i.e., the motion of a single electron in a random magnetic field) by a network with two channels per link and 
local U(2) gauge invariance. Our mapping onto a nonlinear a model easily extends to include this case. The coupling 
constants of the continuum field theory were found to be a x °J = a X y = 2 x 1/2 (mod 1) = at the symmetric point 
of the random flux problem. This is a massive theory with exponential decay of all correlation functions are large 
scales. Thus, contrary to claims made in the literature, there is no room for truly extended states in the random flux 
problem, at least not by slight deformation away from the two-channel network model. This conclusion had already 
been reached in J2^] and Q|. 

Another spinoff helps us along in our quest to understand the integer quantum Hall transition. We observed that 
the SUSY reformulated Chalker-Coddington network model has the structure of what is known as a vertex model in 
statistical physics. Motivated by information from the theory of integrable systems, we then modified the network 
model in several ways. First, we abandoned the alternating direction of the electron's motion along the horizontal and 
vertical straight lines of the network. Instead, we took the direction of motion to be constant along every such line. 
The mapping onto Pruisken's model indicates that this modification changes the physics: the value of the topological 
coupling now is a x °J = 0, which corresponds to a noncritical state. We argued that criticality can be restored by 
doubling the number of channels on all horizontal links (or on all vertical links). Alternatively, we can spatially 
separate the two channels on horizontal links and return to a model with only one channel per link, at the expense of 
doubling the size of the unit cell in the vertical direction. 

The resulting modified network model is critical, but no longer has a justification from a microscopic picture of 
guiding center drift along equipotentials. Its virtue is that the corresponding supersymmctric model is of a type 
for which systematic ways of constructing solutions of the quantum Yang-Baxter equation are in principle available. 
Whether our model can actually be deformed into one that is a Yang-Baxter integrable, is a question whose answer 
lies far beyond the scope of the present paper. 
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